Open charm production at RHIC 
Xin Dong a 

a Department of Modern Physics, University of Science and Technology of China - USTC, 
96 Jinzhai Road, Hefei, Anhui 230026, China 

Recent experimental measurements on open charm production in proton-proton, proton 
(deuteron)-nucleus and nucleus-nucleus collisions at RHIC are reviewed. A comparison 
with theoretical predictions is made. Some unsettled issues in open charm production call 
for precise measurements on directly reconstructed open charm hadrons. 

1. Introduction 

The ongoing four experiments at the Relativistic Heavy Ion Collider (RHIC) are de- 
signed to search for and measure the quark-gluon plasma (QGP), a new state of matter 
composed of deconfmed, locally thermalized quarks and gluons. The equilibrated matter is 
expected to be described by the equation of state (EoS) with partonic degrees of freedom. 
The partonic pressure gradient and the temperature are two important characteristics 
within such kind of EoS. 

The physics results from the first three-year runs at RHIC demonstrate that the par- 
tonic pressure gradient has been developed during the system evolution in heavy ion 
collisions. This has been illustrated in the "white papers" from four experiments [[!]. To 
determine the partonic EoS, the next task is to test the local and early thermalization hy- 
pothesis experimentally. Heavy quark (c, b) is an ideal probe to this end. Due to its much 
heavier mass, it requires more rescatterings to reach the comparable collectivity as light 
quark (u, d, s). If heavy quark collectivity is observed, there must be even more rescatter- 
ings happening among light quarks than expected, because the rescattering cross section 
among light quarks is larger than that between heavy and light quarks. So heavy quark 
collectivity can be used as an indicator for the thermalization of light flavors, although 
heavy quarks themselves do not have to be thermalized [|2|. 

Since charm quark creation requires a large momentum transfer i.e. Q > 3 GeV, it 
is believed to be less affected by non-perturbative effects in pQCD calculations. So the 
measurement on charm quark production in proton-proton (p + p) collisions not only 
provides a necessary reference for heavy ion collisions, but also enables us to test the 
pQCD calculations on both total and differential cross sections. 

In heavy ion collisions, the theoretical calculation shows that charm quarks are mostly 
created through initial gluon-gluon fusions • And since charm quark mass is much larger 
than the estimated system temperature, its production is little influenced by the thermal 
component. Unlike light quark, heavy quark mass is dominated by its current quark mass 
- the mass originating from the coupling with the electroweak Higgs field [H]. Therefore 
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heavy quark is an ideal penetrating probe to the rescatterings and thermalization at the 
early stage of heavy ion collisions. 

The radiative energy loss of charm quark in vacuum is characterized by the "dead-cone" 
effect [E]. When an energetic charm quark traverses through the dense medium, it will 
interact with surrounding partons. Theoretical calculations predict that the suppression 
of the nuclear modification factor (Raa) f° r charm quarks in central nucleus-nucleus (A 
+ A) collisions is smaller than that of light quarks [ HJJ |S] . Most of these predictions 
were made based on the radiative energy loss mechanism and the medium properties to 
our knowledge (gluon density etc.). The interaction between charm quark with medium 
can also be reflected by the charm quark elliptic flow (t^)- The coalescence approach in a 
thermalized medium shows that charm hadrons may obtain a finite v% even if charm quarks 
have a zero t>2 [E|- The charm quark collectivity has been studied in an AMPT transport 
model, and the result shows that a large charm quark interacting cross section is needed 
to produce the magnitude of v 2 comparable with that of light quarks [CHI- Measurements 
of the charm quark collectivity will tell us the degree to which charm quarks interact with 
other partons, and then provide us with pivotal information on the early thermalization 
of light flavors. 

PHENIX and STAR experiments at RHIC have the capability to detect charms. PHENIX 
measures open charms through their semi-leptonic decays: electrons in central arms in 
mid-rapidity and muons in muon arms in forward/backward rapidities. STAR can recon- 
struct open charm mesons through theirs hadronic decays in the TPC. It can also identify 
electrons with the help of other sub-detectors [EJE]- The advantage of PHENIX is that 
it has low budget materials in its inner detectors, while STAR has a large acceptance in 
the TPC around the mid-rapidity. 

2. Charm production in elementary p + p and d + A collisions 

The first reconstruction of open charm hadrons through their hadronic decays was 
reported by the STAR Collaboration in the last Quark Matter conference [ ITT] and was 
recently published in Ref. [E]- The reconstructed charm hadrons and decay channels in d 
+ Au collisions are D° — > K~tt + , D* + — > D°tt + and their charge conjugate channels. The 
Pt coverage is px < 3 GeV/c for D° and 1 < pr/(GeV /c) < 6 for D* + respectively. Event 
mixing technique was used to construct the combinatorial background in the invariant 
mass spectrum. 

STAR also reported the results of non-photonic electrons mostly from heavy flavor 
decays [El El- Results from independent analysis using three different electron identi- 
fication methods (vdE/dx, dE/dx+TOF, dE/dx+EMC) are in good agreement in both 
p + p and d + Au collisions. The non-photonic electron spectrum is also consistent with 
that deduced from STAR's data for open charm mesons. The combined fit of the total 
charm cross section for D° and non-photonic electrons is a^ 5 N = 1.4 ± 0.2 ± 0.4 mb at 
^snn = 200 GeV [E]- Within errors, the electron spectra in p + p and d + Au collisions 
show approximate scaling, implying no significant nuclear effect in d + Au collisions. 

PHENIX reported its spectra of the non-photonic electrons from charm hadron decays 
via three independent methods: the cocktail, the convertor, and the 7 — e correlation [ 
El E|- These methods all give a consistent result: a^ N = 0.92 ± 0.15 ± 0.54 mb in 
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p + p collisions compatible with the STAR data. Measurements of centrality dependence 
of non-photonic electrons in d + Au collisions also shows the N^n scaling [I13|. 
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Figure 1. The total cc cross section per nucleon-nucleon collision vs. the collision energies. 
The low energy data points are selected from fixed target experiments [EH EE]- The 
diamonds depict two cosmic ray measurements [ [T7j . The dashed and dot-dashed lines 
are taken from [IT2"]. 



Available data for the cross production cross section at various energies are shown in 
Fig. [TJ The results in Au + Au collisions at RHIC will be discussed in the next sec- 
tion. The dot-dashed curve depicts a typical next-to-leading (NLO) pQCD calculation 
where parameters are optimized to fit the low energy data The dashed curve is 

a PYTHIA (version 6.152) simulation with the parton distribution function CTEQ5M1. 
Both the NLO pQCD calculation and the PYTHIA prediction give a total cross section 
of 300 — 450 fib, 2 — 3 times lower than the experimental data. Data points from cosmic 
ray experiments also support a large cross section at yfs ~ 300 GeV [Ej. Recent analysis 
of open beauty measurements by CDF shows high order processes {e.g. initial/final radi- 
ation, gluon splitting, and parton shower production) contribute to a large part in heavy 
flavor production at Tevatron [I19j. The discrepancy at RHIC energy between data and 
predictions indicates that these processes may play more important role in theoretical 
model than previous thought. 

Apart from an overall normalization factor, the electron spectral shapes measured by 
PHENIX and STAR are consistent with each other within errors. Shown in Fig. El the 
measured spectra are clearly harder than the overall contribution from charm and bottom 
decays in the PYTHIA simulation. An important issue in obtaining the charm production 
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Figure 2. Non-photonic electron spectrum in p + p collisions compared with PYTHIA 
model LO calculation. The data points are taken from [E] and [El- The PYTHIA 
model parameter setting is inspired by the publication for Au + Au 130 GeV data [ 120] 
with cr cg = 658 ub. 



in the electron approach is how to determine the bottom contamination. A recent NLO 
pQCD investigation at RHIC energy tells us that the crossing point between the electron 
spectra from charm decays and bottom decays may vary in a broad pt range (~ 3 — 10 
GeV/c) [ 121] . This will bring very large uncertainties in electron data from bottom decays. 

The reconstructed D° spectrum covers more than 90% of the total charm yields, while 
the electron spectrum at pr > 0.8 GeV/c only covers ~ 15%. To establish a good 
reference for studying charm production in heavy ion collisions, precise measurements on 
reconstructed charm hadron with large px coverage are necessary. 

3. Open charm production in A + A collisions 
3.1. Charm yields 

Charm yields in heavy ion collisions are expected to be scaled by Nf, in since most 
charm quark pairs are created in the initial hard scatterings. A recent publication from 
the PHENIX Collaboration reported the centrality dependence of non-photonic electron 
spectra in Au + Au collisions at ^/snn = 200 GeV [122]- The electron spectra in all 
centrality bins show approximate Nu n scaling with respect to p+p collisions. The obtained 
charm total cross section in minimum bias Au + Au collisions is 0.622 ± 0.057 ± 0.160 
mb, compatible with the that in p + p collisions. 

At this conference STAR reported reconstructed D° signals in minimum bias Au + Au 
collisions using the same method as that in d + Au collisions [123- They also reported 
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centrality dependence of non-photonic electron measurements in Au + Au collisions [ 
1231 Ej • By combining D° and non- photonic electrons, the extracted total charm cross 
section per nucleon-nucleon collision in minimum bias Au + Au collisions is a^ N = 
1.13 ± 0.09 ± 0.42 mb, consistent with the STAR d + Au result and the PHENIX p + p 
result. The cross section results from Au + Au collisions, measured by PHENIX and 
STAR, are illustrated in Fig. [T] 
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Figure 3. A summary of measurements of the differential cc cross section at mid-rapidity 
per nucleon-nucleon collision vs. the number of binary collisions at ^/snn = 200 GeV. 
The dotted line together with the grey band depict a typical NLO pQCD calculation with 
uncertainties from Ref. [EI]- The dashed and dot-dashed lines are linear fits to all data 
and PHENIX data respectively. 



A summary of available data on charm cross sections at 200 GeV is given in Fig. El 
All measurements agree with each other within large error bars. The dotted line together 
with the grey band display a typical NLO pQCD calculation with uncertainties [EI]- N bin 
scaling is assumed in the extrapolation from p + p to d + Au and to Au + Au collisions. 
The data are systematically above this band. One sees from data that there is a slightly 
decrease from p + p (d + Au ) to peripheral Au + Au , and then to central Au + Au 
collisions, as shown by the dashed and dot-dashed lines which fit to all data and PHENIX 
ones respectively. But we cannot claim much due to large errors. 

The total cross section measurements are important references for charmonium produc- 
tion whose enhancement or suppression in central Au + Au collisions is thought to be a 
robust signal of the QGP. More precise charm measurements in various centralities in Au 
+ Au collisions are needed. 
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One interesting observation, described in Ref. is that the centrality dependence 

of charm yields can be explained marginally by the differential cross section of inclusive 
hadrons integrated over p? > 1-5 GeV/c ~ mp. This means the production of a variety 
of particles is not sensitive to the flavor quantity once the momentum transfer is above 
the threshold. 

3.2. Charm quark energy loss in medium 
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Figure 4. Recent measurements on non-photonic electron Raa i n central Au + Au colli- 
sions from PHENIX (top 0-10%) and STAR (top 0-5% from EMC and top 0-20% using 
d+Au as the reference from TOF) experiments compared with theoretical predictions 
(top 0-10%) from [ 27] (left plot) and [ESI ( ri S ht P lot )- 

Given in Fig.0]are recent results of the nuclear modification factor Raa of non-photonic 
electrons in central Au + Au collisions from PHENIX [HH and STAR [|23J|2I]. The data 
give a consistent, and surprising fact: the suppression factor for non-photonic electrons is 
~ 0.2 — 0.3, which is almost at the same level as that of charged hadrons in the similar px 
range. Two recent pQCD estimations in the radiative energy loss scenario are also shown 
in that figure [ I27J I2H1- These approaches try to fix the transport parameter (dN g /dy 
or q) boundaries by fitting to the Raa for light hadrons The boundaries obtained are 
1000 < dNg/dy < 3500 and 4 < g/(GeV 2 /fm)< 14 respectively. One sees in Fig. H 
the upper limit to which energetic partons lose the largest fraction of their energies in 
the medium due to gluon bremsstrahlung from these two approaches. The comparison 
with the data illustrates the suppression of electrons from charm decays may reach as 
low as that of light hadrons. However, if the bottom contribution is included according 
to pQCD calculations, the overall electron Raa will increase to ~ 0.4 — 0.5. This is a 
significant discrepancy compared to the data at 4 < pT/{GeV/c) < 7. If the data in 
Fig. E] are confirmed to be correct, this will bring at least two open issues: (i) if the 
current radiative energy loss mechanism persists, there is no much room for the bottom's 
contribution in the non-photonic electron spectrum up to px ~ 7 GeV/c. (ii) if the 
bottom's contribution is as what is given by the generic pQCD predictions (the crossing 
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point is ~ 3 — 5 GeV/c), there must be other energy loss effects besides gluon radiation. 
These are challenges to theorists. 

In several recent publications, some authors argued that because in momentum coverage 
jv ~ 1, heavy quark is not ultrarelativistic, elastic collisional energy loss may play an 
important role when charm quarks traverse the medium [EHES]- They computed the Raa 
in the hydrodynamic transport scenario, which gives the strong suppression as observed. 
The approach sheds light on the solution to the present discrepancy. To decouple the 
above two issues, one should precisely measure reconstructed open charm hadrons instead 
of electrons. Certainly, a precise reference in p + p collisions on reconstructed open charm 
hadrons is needed. 




Transverse momentum p T (GeV/c) 



Figure 5. Recent measurements of V2 for non-photonic electron in minimum bias Au + 
Au collisions from PHENIX [ESUS! and STAR [E2j and the comparison with theoretical 
predictions in [ US] 



3.3. Charm quark elliptic flow 

PHENIX and STAR recently measured v 2 for non-photonic electrons in minimum bias 
Au + Au collisions [ESI EH E2]- see Fig. El At pt > 2 GeV/c, even with the claimed 
~ 20 — 30% systematic error by STAR, the measurements are not quite consistent between 
two experiments. In terms of the magnitude, V2 for non-photonic electrons is comparable 
to that of other hadrons [|33| in px < 2 GeV/c. Two model predictions are also shown 
in that figure. Since there is an inconsistency between two experiments and the bottom's 
contribution is uncertain at px > 2 GeV/c, let us focus on the range at px < 2 GeV/c. 
In the coalescence model for a thermalized system [Ej, the picture that charm quark t> 2 is 
the same as light quark v 2 is supported by the data. Compared with the transport model 
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results [llOj. the data favors charm quark has a large rescattering cross section, indicating 
that the charm quark has a finite Although it is hard to extract v 2 of charm hadrons 
from electron t> 2 , the data suggest there can be a non-zero v 2 for the charm quark. If 
the non-zero charm quark collectivity is confirmed, it means light flavor thermalization 
at RHIC, as I have argued in the introduction section. 
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Figure 6. Raa of charm hadrons (thin lines) and decayed electrons (thick lines) for several 
assumed D spectra in central Au + Au collisions. The charm hadron spectrum reference 
used in p + p collisions is from STAR's d + Au measurements [EJ. The thin dashed line 
depicts the Raa of charm hadrons same as that of charged hadrons. The thin dot-dashed 
line is for charm hadrons with the Blast-Wave (BW) behavior at Tf = 160 MeV and 
(/3t) = 0.4c in low pt and with the same behavior as charged hadrons' Raa in high pt [ 
P. The thin solid line is for D with the same BW parameters as the dot-dashed line in 
low p T , but the total yield of charm hadrons is assumed to obey N bin scaling. While in 
high pt, their Raa is taken from the dead-cone energy loss calculation in Ref. [E]. Thick 
lines depict the Raa of electrons from charm decays with the distribution of the same line 
style. 



3.4. Complementary remarks 

From the measurements of Raa and v 2 for other particles, the radiative energy loss 
can contribute to part of v 2 , but not much. In the recent work mentioned above [ 29J, 
elastic collisions may provide a large fraction of energy loss of heavy quarks when ~ 1. 
During hydrodynamical evolution, the flow of underlying medium will influence the heavy 
quark spectrum and heavy quark will pick up some flow. In this case, Raa and t> 2 are 
quite correlated, which means charm quark Raa must be strong suppressed if a large v 2 
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is observed. So the combination of measurements on the charm quark spectrum and v% 
is essential, especially in low px region. The above arguments are based on the hydro 
assumptions. 

To test the medium response to heavy quarks, the low pt part is quite relevant. How- 
ever, the electron spectrum from charm decays cannot disentangle different shapes in this 
Pt region due to smearing of the decay kinematics [Ell- This effect can also be reflected 
on the Raa of non-photonic electrons [HEj. The simulation results of various input charm 
hadron's Raa (thin curves) and the corresponding Raa for electrons from charm decays 
(thick curves) are shown in Fig. |H1 One sees that although at high p?, the electron Raa 
can reflect the suppression of charm hadrons (solid curves and dashed curves), the elec- 
tron Raa at low pt cannot tell different thermal shapes (dashed curves and dot-dashed 
curves) : the thick dashed and thick dot-dashed curves are almost identical. Therefore we 
need in the future precise measurements of the spectrum for reconstructed open charm 
hadrons at low p? ■ 

4. Conclusion and outlook 

The heavy flavor program has started at RHIC extensively. Heavy flavor collectivity 
is expected to be an ideal probe to the light flavor thermalization. Plenty of new and 
surprising results on open charm production have been presented on this conference. 
However, most of present measurements use electrons decayed from charm hadrons which 
brings large uncertainties. The electron approach can be only a placeholder. We need 
precise measurements for reconstructed charm hadron spectra and v-i in a wide p? range. 
So the current sub-detector upgrade proposals in pipe for PHENIX and STAR detectors 
are very important to this goal [EI3EIEEI1- With the upgraded inner tracking detector, 
PHENIX and STAR can reconstruct the secondary vertices of open charm decays with 
much lower background. We look forward to more exciting physical results from RHIC 
with upgraded detectors in the future. 
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